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SATELLITE-TO-SATELLITE TRACKING 

F.O. Vonbim 
P.D, Argeriticro 
P.E. Schmid 

ABSTRACT 

Tltc possibility of iitili/dng geostationary spacecraft as orbiting tracking 
stations was first considered by NASA’s Goddard Space Flight Center in 
the early 1 960’s, However, satellitc-to-satcllite tracking did not become 
a reality until 19 April 1975 when the geostationary Applications Tech- 
nology Satellite-6 (ATS-6) began to relay Doppler data from the near- 
Earth geodynamics satellite, GEOS-3, Shortly thereafter ATS-6 radio 
tracking of the NlMBUS-6 weather satellite was initiated. This paper 
presents results from these satcllite-to-satellite experiments in terms of 
tracking system performance and orbit determination accuracy. The ex- 
perience gained during these tests is directly applicable to the NASA 
geostationary Tracking and Data Relay Satellite System (TDRSS) which 
has a first launch scheduled for mid 1 980. 
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SATELLITE-TO'SATELLITE TRACK ING 
1.0 INTRODUCTION 

The possibility of using geostationary satellites for communications was discussed in 
the popular literature as early as 1956 (1), The first detailed jjroposal for a synchronous 
tracking satellite system for the purposes of orbit determination was provided by Von Bun 
in 1967 (2, 3). Since then a number of papers (4, 5, 6, 7, 8, 9,) have considered the use of 
satellitc-to-sateilite tracking for orbit determination and for gravity field model refinement. 
These papers mention that with regard to coverage, a satellite-to-satellite tracking system 
has a significant advantage over ground bused tracking systems. For instance, with a singic 
synchronous relay satellite, a satellite-to-satellite tracking system is capable of observing an 
earth orbiting satellite during almost half of every orbit. Equivalent coverage of a satellite 
in a high inclination orbit would be difficult to obtain with a ground based system. 

In 1 968 during the early planning phases of the geostationary ATS-6 and near-Earth 
NIMBUS-5 experiments it became clear that this satellite configuration would be ideally 
suited to evaluate the concept of satellite-to-satellite tracking and to provide valuable ex- 
perience in processing this new data type. The experiment as defined in October 1 968 ( I Ol 
incorporated both radio time delay (range) and Doppler frequency shift (range rate) mea- 
surements. This experiment, entitled the “Tracking and Data Relay Experiment” (T&DRE) 
was conducted as planned except that NIMBUS-6, which was launched June 12, 1975, 
rather than NlMBUS-5 carried the T&DRE equipment. In early 1972 plans were completed 
for a very similar ATS-6/GEOS-3 satellite-to-satellite tracking experiment. The GEOS-3 
satellite was launched on April 9, 1975. Another satellite-tosatellite tracking effort involv- 
ing the ATS-6 was the Goddard Apollo-Soyuz Geodynamics Experiment (1 1) performed 
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during 1 975. However the accent of this experiment wiis gravity anomaly detection rather 
tlian orbit determination, Tiic ATS-6, wlileli was tlic relay satellite for tiicsc experiments, 
was laundied on May 30, 1974 and is still in operation. 

The results of these experiments are relevant because NASA intends to use tlic Track- 
ing and Data Relay Satellite system (TDRSS) (12) for operational orbit determination of 
NASA satellites. The system will consist ot two synclironous relay satellites (one at 41 
degrees west and one at 171 degrees west) and a common ground station under construction 
at White Sands, New Mexico. Operations will begin in November 1980. Hence by the early 
nineteen eighties satellitc-to-satellitc tracking data will be routinely processed to obtain 
orbits. 

This paper is a report on the results of the ATS-6/GBOS-3 and the ATS-6/NIMBUS-6 
satellitc-to-satellite t raeking orbit determination experiments. The tracking systems used in 
these experiments differ from the TDRSS, primarily in the use of one rather than two syn- 
chronous relay satellites. However tlic authors believe and simulations mentioned in this 
paper indicate that the insights gained from the experiments with regard to proper data re- 
duction techniques and expected results arc applicable to ilie TDRSS. 

1.1 EXPERIMENT SPACECRAFT 

Tlie key to all satellitc-to-satcllitc experiments to date has been the geostatlonaiy 
ATS-6 spacecraft (13, 14). During the past three years the equatorial ATS-6 has been sta- 
tioned over both tlic Pacific in proximity of continental U.S.A, and Africa. Accordingly, 
ATS-6 ground stations have at various times been operated at Rosman, North Carolina; 
Mojave, California and Madrid, Spain. The ncar-Earth satellites tracked via ATS-6 have been 
GEOS-3 ( 15), Apollo-Soyu7.( 16). and NIMBUS-6 (17). 

The nominal GEOS-3 orbital parameters arc a mean altitude of 843 km, an inclina- 
tion of 115°, an eccentricity of 0.004, and a period of 101 .8 min. The orbit parameters 
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were chosen to minimize resonance of the subsatcllite trace with any given Earth feature 
and to provide orbit traces which cover the Earth in a gridwork pattern. 

Tlie ApollO'Soyuz mission included the Geodynamics experiment where Apollo was 
tracked via ATS-6 for the mission duration (15 July to 24 July 1975). The nominal Apollo 
orbit at insertion was 150 km by 170 km at an inclination of 51.8^. 

Finally, the NIMBUS-6 weather satellite is in a Sun synchronous polar orbit with a 
mean altitude of 1 1 10 km, an inclination of 100°, and a period of 107.4 min. 

2.0 SATELLITE-TO-SATELLITE TRACKING 

A satellite radio or laser tracking system makes measurements of such parameters as 
range, range rate, angles and direction cosines to a spacecraft relative to a given tracking sta- 
tion. In two-way tracking a signal is transmitted from a well surveyed ground station to a 
spacecraft transponder wiilch frequency translates the signal for re-transmission directly 
back to the ground station or, as in the case of satcllitc-to-satcllite tracking, to another 
spacecraft. The two-way tracking system developed for the experiments discussed in this 
paper measures "range” in terms of the round-trip time delay on a 100 kHz tone and range 
rate in terms of the Doppler shift on a 2 GHz carrier signal (14, 18). 

2.1 GEOMETRY 

The tracking geometiy is shown in figure !. The ground station transmits a signal to 
the near Earth satellite via the synchronous spacecraft. This same signal is “turned around” 
and transmitted (at a slightly offset frequency) back to the ground site again via the high 
altitude satellite. For purposes of stability NASA geostationary orbits have been maintained 
at inclinations which extend from 1,5° to 6°. As a consequence the path indicated as Rj 
(figure 1) varies as a function of time as docs R 2 (13). Because of the radio propagation 
times Involved and the fact that both spacecraft are in motion relative to the ground site, 
four distinct paths must be considered when interpreting the Doppler (range-rate) and time 
delay rangc/measurernents (13, 19). 
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riie “range’* iivasurenient is pcrl»)rmcil by cuinpuring Iransiiiitleii and received t«>ne 
zero crossings, the Highest resolution tone fret|uency in this case being 100 Kll/. Lower Ire- 
i|uency tones are SL't|uentially used during actpiisition for ambiguity resolution. Ihe lower 
tones are at 20 Kll/. 4 KH/. KOO II/.. loO II/.. .12 II/. and H II/. 

I'he tone ranging measurement is(|uite straightldrward and ranging accuracy depends 
chiclly on the »|uality of preHight calibration of both the A I S ami NIMHL’S transponder 
group delay. Such pretlight calibration data have Iven taken over a range of fre(|uencies and 
temperatures. Indications are that with careful calibration the total systematic delay error 
in the ranging measurement can be held to a few meters of c(|uivalent one-way range. 

I he “range rate” measurement is performed by counting cycles of Doppler over a 
measiircil time interval In the case of Cd OS-.l and .Apollo .’he measurement consisted of 
the number of Doppler cycles accumulated m the regular sampling interval ( I or 10 seconds 
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dcpcmliin; on mode of operation). For NIMBUS-6 tlie incasiircmcnt consisted of the time 
interval required to accumulate a fixed number of Doppler cycles ( 1 8), The electronics 
for ATS-6 satellltc-to-satcllitc tracking have been so configured that tlic Doppler output is 
approximated by: 

+^2)) 

where 

fj = jneasured average Doppler frequency 
fj = uplink frequency 
c = ,<:peed of iiglit 

k, a| and are scalar constants determined by equipment frequency multiplications 
f j = average range-rate ATS to ground site 
T 2 = average range-rate ATS to NIMBUS, Apollo, or GEOS-3 
A detailed discussion of Doppler factors in satcllitc-to-satcllite tracking is given in 09). 

Tiic upiink to ATS-6 (f^) is at a nominal 6 GHz, The link to and from the low satellite is 
nomiiially 2 GHz and ATS-6 back to ground at 4 GHz. The range and Doppler measure- 
ments will als<: be biased by the Earth’s troposphere and ionosplicre. Measurement biases 
up to meters in range and tens of cm/scc in range rate can be expected at 2 GHz. Atmos- 
phere refraction effects can to a large extent be modeled out. Some of the work done in 
this area at NASA-GSFC is indicated in (20, 21 , 22), The atmospheric range bias is fre- 
quency independent through the troposphere and inversely proportional to frequency 
squared through tlie ionosphere. The range rate bias, in addition to the foregoing, is pro- 
portional to the rate of scan through the atmosphere as well as to the magnitude of hori- 
zontal gradients, 
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2.3 ORBIT DliTBRMINATION TECHNIQUES 


Problems of Orbit Determination witli SateIlitc*to*Satcllltc Tracking 

The unfamiliar feature of iletcrmining user satellite orbits by means of a satelliie-to* 
satpllile tracking system is the presence of the relay satellite state as an error source. The 
simplest procedure for estimating user satellite state in the presence of this error source is to 
estimate a satellite epoch state from the satcllite-to-satcllltc tracking data with the relay 
satellite state constrained to a irreviously determined orbit and left unadjusted in the rcduc* 
lion process. With this approach the uncertainty in relay satellite state is manifested as an 
uninodcicd and time varying error source which alters the estimate of user satellite state. 
Some subtleties arc encountered in determining the effect of this error source. The time 
history of relay satellite state error is a function uf the way in which the epoch state was 
computed, Eor example, suppose the relay satellite is independently anti continuously 
tracked over a given period and a least scpiarcs algoritlmi used to estimate epoch state at the 
beginning of the period. If this epoch state is then propagated to the enti of the period 
using the same dynamic model that was used to jrrocess the data, the resultant errors will be 
constrained by the data fitting criterion implicit in tlic least squares reduction algorithm. 
Tile errors so obtained will be smaller than the errors obtained if either one did not match 
tiynamic models or if one propagated the epoch state beyond the data collection period, 

The same phenomenon can be understood from a statistical vantage point by observing that 
when the dynamic models arc matched the epoch slate , errors become correlated with dy; 
namic parameter errors, and that over the data arc these correlations tend to minimize the 
errors in the epoch state propagation. 

Relay satellite stale uncertainty appears to be a significant error source even when the 
relay satellite or satellites arc continuously and independently tracked. Argenticro and 
Loveless (23) simulated the orbit recovery of a satellite in a 300 km, polar, circular orbit 
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with the Tracking Data Relay Satellite System (TDRSS) (24). The TDRSS satellites can 
relay range and Doppler infomiiatlcn from a low altitude user satellite to a ground station. 
Tlie simulations assumed that each synchronous satellite was continuously tracked from 
two ground stations and that 24 hour data spans were processed to estimate user satellite 
state. The same dynamic models which were employed to estimate relay satellite epoch 
states were also used to estimate user satellite state from the satellitc-to*sateliitc tracking 
data. The effect of Geopotential and atmospheric drag errors were included in the simula- 
tion. The results showed that user satellite position could be recovered with an average 
total position error of 260 m. Tlic major part of tliis error is caused by the error in esti- 
mates of relay satellite cpoclj states. When these simulations arc repeated without the as- 
sumption of continuous tracking tlie results arc considerably worse. 

A standard approacli to dealing with troublesome error sources in an orbit determina- 
tion is to augment tlie list of estimated parameters in the data reduction by including these 
error sources. This approacli can certainly be implemented with regard to relay satellite 
state errors by simultaneously estimating user and relay satellite epoch states from informa- 
tion supplied by the satcllite-to-satcilitc tracking data. From one vaniage point this is an 
undesirable solution in that the user is uninterested In the state of the relay satellite and 
would rather not burden the numerical procedures witli tlie need for simultaneously esti- 
mating relay satellite state with the user satellite state, However, the results of independent 
covariance analyses performed by Fang and Gibbs (25), and Argentiero and Garza-Robles 
(26) indicate that an unconstrained simultaneous estimate of user and relay satellite states 
using satcllitc-to-satcllite tracking data can yield an estimate of user satellite state which is 
consistently better than 100 m. 
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Numerous sluiultflucous unconstrained solutions have been attempted using range 
sum aiul range sum rate measurements obtained from tlic ATS-6/CJliOS*3 combination and 
the ATS-6/NIMUUS'6 combination atui in all case# the solutions have been Inaccurate and 
numerically unstable, Clearly our c.'tperlcncc with real reductions of salcllite*to*satclllte 
Hacking data is not compatible with the results of previous error studies. In order to under- 
stand tlic discrepancy we have performed a junncrical simulation of the ATS-6/GliOS-3 
satcllite-to-satellite tracking experiment. The difference betwccjl a numerical simulation 
ami a covariance analysis can be described as follows: in a simulation, data are generated 
and a least squares adjustment process is actually performed. The estimated state is then 
comimrcd to the reference or unperturbed state at various points along the orbit and conclu- 
sions can be drawn concerning the accuracy of the process. In a covariance nnalysis mode, 
the least squares adjustment process is jmstulated rather than actually performed, and under 
the assumption that 'v.v*'* ‘he range of expected errors, perturbations of orbital estimates arc 
approximately linear functions of perturbations of the error .sources, the associated covari- 
ance matrix of the epoch state recovery Is computed. Willi the aid of state transition matri- 
ces the covariance matrix at epoch can be propagated to obtain the covariance matrix of the 
satellite state recovery at any point in tlie orbit. 

For tlie numerical simulation a computer program was used to generate 1 2 hours of 
range and Doppler satcllitc-to-satellile tracking ciata from Hie ATS-6/GEOS-3 satellite com- 
binaiion. In tliis data generation the Naval Weapons Laboratory (NWL) geopotential field 
was used, A random number generator added wliitc noise of standard deviation 1 mm/sec 
to the Doppler data and wliite noise of 2 m to tlie ranging data, values consistent with track- 
ing system performance. The SAG 69 gcopotcntial field and an orbit detennination pro- 
gram were used to reduce tlie data to simultaneously estimate the ATS-6 and GEOS-3 cpocli 
states, 'file estimates GEOS-3 epocli state was propagated along tlic entire 12 hour data 
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collection period using the SAO 69 gcopotcJitial field. This orbit was compared at selected 
time points to the true CnOS'3 orbit which was obtained by propagating theGEOS-3 ref- 
erence epoch state with the NWL geopotentiai field. The average difference between the 
two orbits was over 900 m. Also the nominal covariance matrix of the data reduction re- 
vealed that several correlations between estimated parameters were of absolute value near 
unity. This implies that the normal matrix which is inverted in the least squares estimation 
process Is poorly comlitior.ed. l ienee small perturbations of the elements of this matrix 
such as those caused by computer roundoff and other effects cause major perturbations 
of the elements in the inverted matrix. This amplification effect in the inversion of a poorly 
conditioned matrix can lead to an inaccurate estimate of a satellite epoch state or in some 
eases a divergence of the least squares interatlon procedure, This is the probable cause of 
poor results using a simultaneous esvimation approach in both the simulated and real data 
reductions. In a covariance analysis of the simultaneous estimation approach the least 
squares algorithm is not actually executetl and consequently these numerical problems arc 
never manifested. For this reason the techniques of covariance analysis provide a somewhat 
optimistic assessment of orbital accuracies obtainable from simultaneous estimation with 
satcllite-to-'^atcllltc tracking data. 

Thus the two conclusions of our analyses arei 1) The uncertainty in relay satellite 
state is a significant error source which cannot be ignored in the reduction of satellite-to- 
satellitc tracking data and 2) that based on both simulations and real data reductions it is 
nunierically impractical to use simultaneous unconstrained solutions to determine both 
relay satellite and and user satellite epoch states. The estimation technique used to generate 
the results shown in subsequent sections may be described as a Bayesian or least squares 
with a-prlori procedure. This approach permits the adjustment of relay satellite epoch state 
in the reduction of satellite-to-satellite tracking data but without the numerical difficulties 
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introiliiccc) by an ilNconclitioncd nonnni matrix. Theoretically this technique obtains the 
best possible estimate of user satellite state based on all available information. A mathemati' 
cal description follows. 

Mathematical Description 

In this mathematical development we assume the existence of two separate data sets: 
yj - Ranging observations between ATS-6 and ground based tracking stations 
V 2 - Satellitc-to-satellite tracking of user Sittellitc (range sum and range rate sum) with 
ATS-6 as relay satellite. 

The parameter set to be estimated consists of two satellite epoch states. 

X j ~ Six dimensional ATS-6 state at epoch time Tj 
X 2 ” Six dimensional user satellite state at epoch time T 2 

Tii.p data set yj is corrupted by errors in the measuring process, Hence represent yj as: 

yj “V] +Vi,c(Vj)='0,c(vjV,'^) = Q, (1) 

where yj is the correct or noiseless representation of the data set, vj is a vector of random 
errors of zero expectation and covariance matrix Qi . Describe the functional relationship 
between and X| as 

yi=f(X|) (2) 

The right .side of eq. 2 represents a computational algorithm obtained by integrating satel- 
lite motion to each observation time and coniputing the ideal observations, The standard 
Ir'jr.t squares estimator x j of x j is that vector which minimizes the loss function, 

Uxj) = (y, -f(x,))'>’Q,-'(yi-f(x,)) (3) 

Assumiirg the linearity of eq, 2, the vector which minimizes the right side of eq, 3 is also 
known to be a minimum variance estimator. A first approximation to the desired minimum 
can be obtained by expanding equation 2 in a first order Taylor series about nominal value 
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whcr'^ 6 y, ami fix, arc deviations of and x, from nominal values and A, is the so-called 
sensitivity matrix. The estimate of fix, is 

fix, =(A,TQj-lA,rU,TQj“l6y, (5) 


where 


5 yi=yi“f(x,,„) 

The vector fix, is added to x , „ to form an estimate of x, . This estimate can be used as 
a new nominal and the process can be repented until a convergence criterion is satisfled. 
The covariance matrix of the least squares estimate x, of x, is 


c = eax, -X,] [xj -x, ]'*') = (aT’Qj-I A)-* (6j 

The next step is to obtain an optimal processing of the data set y2- Define a 12 dimensional 
vector z as 


z = 



( 7 ) 


Represent the data set y2 as 


y. ••= y2 + v,, efvj) = 0, e(v2V2'^>’-- Q 2 


( 8 ) 


where y2 is the correct or noiseless representation of the data set, V2 is a vector of random 
errors of zero expectation and covariance matrix Q2. The functional representation be- 
tween y2 ^nd zis presented as. 


y 2 = e(z) 


( 9 ) 
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As was tiic ciise with equation 2, the right side of eq, 9 represents a computational algorithm 
involving the integration of satellite equations of motion, 

Tile least squares estimate of z would not be optimal unless all available information were 
included in the loss function. Hence it is appropriate to treat the least squares or minimum 
variance estimate x j of x j as an a-priori estimate weighted by the inverse of the covariance 
matrix provided by equation 6, The resulting loss function to be minimized is 


Uz) = (y2-g(z)rQ2"'(y2-8(2'0) 


( 10 ) 


Again, the required minimum can be obtained iteratively by expanding equation eq. 9 in a 
first order Taylor scries about a nominal estimate of Z 


ASA Qe(z) 
5y2 - A25 z, A 2 - 


Z-Zr 


01 ) 


where ^2 deviations of y 2 and z from nominal values and A 2 is the sensitivity 

matrix. The estimate of Sz is 


whui'e 




5xj = Xj - Xj , fiy2 = yj - g(zii) 


( 12 ) 


The vector 8z is added to Zn to estimate z. This estimate is used as a new nominal and the 
process is repeated until a convergence criterion is satisfied, The final covariance matrix for 
the estimate of satellite state x j and satellite state X 2 is 

/ \-l 


E 


r/. 



'^ 1 ' 

VI 

Lr 


)r 


lj J 


- ( A2^Q2“^A2 + c"1 


(13) 
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It can be shown that tlic two step process defined above in which data set yj is processed 
and then data set ,V 2 is processed is equivalent to a single step unconstrained least squares 
estimation of X] and X 2 using both data sets y | and y 2 > Hence this procedure leads to the 
most accurate estimate of both user and relay satellite state based on available information. 

3.0 EXPERIMENTAL RESULTS 

The experimental results can be considered in three categories -- namely, 

• tracking system performance 

• geostationary satellite orbit evaluation 

• near Earth satellite orbit evaluation 

3,1 TRACKING SYSTEM PERFORMANCE 

The expected error for the NASA range and range rate satellitc-to-satellite tracking 
system is a function of many controlled parameters such as range tone frequency, sample 
rate, bandwidth settings, signal-to-noisc spectral density ratios, spacecraft dynamics, etc, 
(13). However, the system is generally used with what might be termed a standard set of 
options such as: 100 kHz maximum range tone frequency, signal levels such that system is 
not thermal noise limited, 1 per second or 6 per minute data rate, and a 25 Hz range track- 
ing loop two-sided noise bandwidth, Tabic I lists the theoretical system performance for 
the foregoing selected options. Doppler averaging time is approximately one half the sample 
time interval for NIMBUS tracking and equal to the sample interval for Apollo and GEOS 
tracking. 

For averaging times, T, up to about 10 seconds the noise decreases as 1/T. The prin- 
cipal Doppler noise contribution comes from receiver voltage controlled crystal oscillators 
and the analog to digital conversion. For longer integration times the Doppler noise is also 
influenced by noise falling off as l/T, an effect attributed to the phase jitter in the 
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Tracking System Measurement Resolution 


Range (Meters) 

Range Rate (Cm/Sec) 

Systematic 

Random 

Systematic 

Random 

1.2 

1.2 

Negiigiblc 

0.03 



NOTE: 10 Sec. Averaging 


transmitter reference signal used at the Doppler extractor. For satellite-to-satellitc tracking 
involving ATS-6 the range-rate resolution for T seconds of iivcraging (27) is given by: 




(t) ' i^)' 


This range-rate resolution versus Doppler measurement averaging time is plotted in Figure 2. 

It should be mentioned that the least significant range bit recorded is 1.5 meters 
which is consistent with the best expected one way performance of 1.7 meters resolution. 
Measured results indicate close agreement with expected system performance. System 
random errors or “noise” arc generally observed by the least squares fitting of short data 
spans (i.e., I to 10 minutes) with polynomials of at least 5th degree to account for space- 
craft dynamics. Care must be taken such that the polynomial itself does not introduce 
apparent error. If tlie data is from a static or collimation tower test a least stpiares straight 
line fit is appropriate. 

Assuming reasonable tracking geometry the accuracy of spacecraft position and 
velocity determination will be primarily limited by tracking system performance for any 
computation spanning the data collection interval. That is^ if continuous tracking is pro- 
vided from a set of well surveyed stations the computation is essentially one of geometry. 
On the other hand the accuracy of orbit prediction based on an initial spacecraft vector 
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Figure 2. Range Rate Resolution 


determination will be degraded as a function of time in direct relation to the accuracy to 
which physical parameters are modeled. This modeling includes gravitational fields, atmos- 
pheric drag and refraction effects, solar pressure, station location determination and so on. 
The most critical of these modeling parameters in terms of orbit determination accuracy is 
the gravity field model which at present is generally expressed in terms of a spherical har- 
monic expansion. 

3.2 ATS-6 ORBIT DETFRMINATION RESULTS 

The optimal reduction of satellite-to-satellite tracking data to determine a user satel- 
lite orbit requires an accurate a priori estimate of relay satellite state, lienee, it is important 
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to determine the expceted errors in tlic estimate of ATS-6 state from ground based tracking, 
More prccisciy, we require answers to tiicsc questions: 

1) How accuraleiy can tlie ATS-6 orbit be determined over an orbital period 
(24 lir) from data wliiclv spans tlie orbital period? 

2) Once an ATS-6 epoci) state is determined, how accurately can that state be 
propagated beyond tlie data arc wliicli was used in its estimation? 

GEOSTATIONARY SATELLITE SHORT 1’ERM ACCURACY 

The first question was investigated by examining reductions of ATS-6 trilateralion 
tracking obtained on November 3, 1974. Trilateration riata is obtained by sending a signal 
from a single tracking station to several strategically deployed unmanned low cost trans- 
ponders via the satellite wliose stale is to be determined. Tlie time required for the radio 
signals to complete the round trip to and from eacli transponder is measured at tlie trans- 
mitter site. Tile interrogating sites were located at Rostnan, Nortli Carolina and Mojave, Cali- 
fornia. The transponders were located at Rosnian, Mojave, Grecnbelt, Maryland, and San- 
tiago, Chile. 

The metliod of “orbit overlaps” was used to evaluate the orbit determination accu- 
racy of tile system. Tin's procedure can be outlined as follows: 

1 ) Determine a satellite epoch state using cacii of two independent data sets 

2) Propagate estimated epoch states over a common or overlapping interval 

3) Difference the two orbits over tlie common interval (differences are usually 
displayed in along track, cross track, and radial components). 

in some cases the orbit overlap metliod can lead to an under-estimation of orbit errors 
since biases in orbit estimates may cancel in orbit differences, Hence, the metliod sliould be 
viewed as a test of the internal consistency of an orbit determination process ratlier tlian an 
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TRANSMITTER TRANSPONDER 


absolute measure of accuracy. Data set I used in tlie orbit overlap test was obtained with 
Rosman as a transmitting site and with transponders at Rosman, Mojave, Grecnbelt, and 
Santiago. Data set 2 was obtained with the same transponder sites but with the transmitter 
located at Mojave. The tracking schedule is shown in figure 3. The two interrogating sites 
are identified in figure 3 under TRANSMITTER as Rosman, North Carolina and the Mojave, 
California "Hybrid Transportable" station. Each data stretch was approximately 5 minutes 
iong and the data rate was one sample per 10 seconds. Separate orbit arcs were computed 
from data set 1 and data set 2. The total position differences between the two orbits over 
the 24 hours of Nov. 3, 1974 were computed and arc displayed in figure 4. The mean posi- 
tion error is about 100 m. A typical .set of range residuals is shown in figure S. The range 
residuals over this arc arc on the order of 20 m. 
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Figure 3. Tracking Schedule 3 November 1974 


Assuming that there are no significant biases in the trilateration orbit determination 
whose effects cancel in the orbit overlap test, the results of figure 4 suggest that continuous 
tracking of ATS-6 over a 24 hour period leads to an orbit estimate over the period which is 
accurate to about 100 m. 
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F'igure 5. ATS-<> Orbit Kosidiiuls 


GEOSTATIONARY SATELLITE LONG-TERM ACCURACY 


lit gcncnil one cunnot ussimie tliat relay satcliites arc continuously tracked. Hence, in 
tiic reduction of satellite>to-satelIite tracking data, it may be necessary to use an estimate of 
user satellite state obtained tliroiigli a propagation titat was unconstrained by tlic data fitting 
criterion of a least squares uigoritlim. Wlicn this occurs the accuracy of the orbit estimate Is 
entirely dependent on the correctness of the force models used in the propagation. 

The orbit overlap technique (30), utilizing data obtained during July 1975 was used 
to estimate the accuracy of a free or unconstrained propagation of an ATS-6 epoch state, 
The data sets used in the overlap tests were*. 

Data Set I — 24 hours of data over July 13, 14, 1975, Tracking stations located at 
Madrid, Ascension Island, and Johannesburg. 

Data Set 2 — 24 hours of ranging data over July 25, 1975. Tracking stations located 
at Madrid, Ascension Island, and Johannesburg. 

Each data set was processed to estimate an ATS-(5 state vector for epoch time July 
16, 1 975 at 7 hr., 25 min. The epoch states were propagated forward for 10 days and along 
track, cross track, and radial differences were computed at 1 5 minute intervals. The root 
mean square along track difference was over 2 km. Figure 6 is a plot of these along track 
differences. 

The large errors which occur during the free propagation of an ATS-6 epoch vector 
must be caused by a misrepresentation of force models. The obvious candidates arc: 

I) Unmodeled venting and thrusting of ATS-6 to accomplish satellite attitude cor- 
rections. Motions due to antenna maneuvering may also introduce errors. 
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2) Mismoilclin(( oT solar ruiiiation pressure. In all data reductions, ATS-T» was as- 
sinned to present a constant cross section to the sun. In fact, this is not the ease. 

3) An error in representation of the central force term of the Tarth’s gravity field. 
An estimate o*' l!ie uncertainty in estimates of this parameter is one part in l(/*. 

Frror source number appeared to us as the most likely cause for the major part of 
the errors exhibited in figure 6. In order to measure the effect of uncertainty in the gravity 
field parameter on the free pro|nigation of ATS*f>, the following simulation w;*s performed; 
ranging observations to the ATS-(i from sites at Kosman, Santiago, and Mojave were gener- 
ated for a three day span. The obsenations were corrupted with white noise with a standard 
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deviation of iO tn. Tlic value of tlio gravity ficid parameter used to generate tlic data was 
perturbed by one part in iO^ and this value was used along with a least squares estimator to 
estimate an epoch state at the beginning of tiie three day data span. Tlie perturbed value of 
tlic gravity field parameter was used to propagate tills epoch state for six days, Over tlie 
three days covered by data tlie propagated orbit differed from tiic assumed true orbit by 
about 100 m. But at tiic end of tlie six day propagation period tlie errors were approxN 
niatcly 2 km, Tlie resuits of tills shiniiation suggest tliat tlie error in tiic central force term 
of tlie Bartli’s gravity field is sufficient to account for tlie errors in the ATS*6 free propaga- 
tion as iminifested in figure 6. 

SUMMARY OF RESULTS 

Overlap tests performed with real data togetiier witli simulation results suggest tliat 
by ])rocessing data over one ATS-6 orbital period, tlie ATSt6 state over the orbital period 
can be detcniilncd witii an average accuracy of about 100 ni. But other results show that 
when longer data arcs are used or when an estimated ATS-6 epocli state is propagated well 
beyond the data arc used in its estimation, errors in tiic kilometer region arc encountered. 
These facts indicate that there arc significant errors in tlie models of the forces acting on the 
ATS-6. The most likely candidate is the cnor in representation of the central force term of 
the gravity field. 

3.3 GEOS-3 ORBIT DETERMINATION RESULTS 

The GEOS-3 orbit determination results were derived from data obtained over the 
weekend of May 3, 1975. The tracking sclicdulcs and the trucking systems used in the eval- 
uation arc sliown in figure 7. The figure shows that five passes of range sum and range sum 
rate data were available. A Bayesian estimation technique described in a previous section 
Was used to obtain two separate and overlapping GEOS-3 orbits. A GEOS-3 epoch state at 
May 2, 22 hr was estimated using iill tlie ATS-6 ranging data and the first three passes of 
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raii)!c sum aiul ran*.!C sum rale ilatu. I he A I S-li raiitimi: data was wei|!liled aecording to a 
standard deviation of 2 m. ( he raii|ie sum and range suir rate data were weighted according 
to standard deviations respectively of 2 m and I mm/sec. I he complete (II.M*7 geopolen* 
tial field was used in this and all other data reductions. Hie estimated epoch state was 
propagated to the end ot the data span of May 22 lir. I he process was repeated with the 
last lour passes of range sum and range sum rale data to estimate a (il OS-3 epoch stale at 
May 3, lU hr. I his epoch state was propagated to the end of its data span at May 4, 10 hr. 
(he total position difference between the two orhils during the 1 2 lir overlap pericnl as 
shown in figure H varies periodically luMween 10 and 25 meters. 

As mentioned in a previous section, orbit overlap procedures can provide an overly 
optimistic asH’ssment of orbit determination accuracy. A more objective measure of accu- 
racy is obtained by comparison with an t)rbil derived from an iiulependent and well 
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calibrated data set. Figure 7 displays tlic C*band tracking available from Wallops Island and 
Bermuda during the weekend of May 3, 1975. A three-day GEOS-3 arc was derived from 
the C-band data and comjjarcd to u similar arc derived from the five passes of satellltc-to- 
satellite tracking data and ATS-6 rangiiig data. The root mean sijuare differences in tlic two 
arcs were: 

radial 5 m 

cross track 200 m 

along track 39 m 

Various orbit results indicate that total position error for C-band derived GEOS-C orbits is 
on tlic order of 50 m (31). Hence, it is only in the cross track direction that the orbit deter- 
mination derived from satellitc-to-satellitc tracking data differs significantly from the 
C-band orbit. The large cross track .'rrors can bo explained in terms of the tracking geome- 
try. For each of the five satellite-to-satellitc tracking passes shown on figure 7, the GEOS-3 
satellite passed almost directly under the ATS-6 satellite. Consequently there was little 
cross track information in tlie range sum and range sum rate data. It is a reasonable assump- 
tion that with a better geometric distribution of passes the cross track errors would be sub- 
stantially reduced. 

3.4 NIMBUS-6 ORBIT DETERMINATION RESULTS 

The NIMBUS-6 overlap results were derived from data obtained over the weekend of 
Feb 8, 1976. For tliis experiment a highly accurate reference orbit suitable for the purpose 
of comparison was unavailable. This Imph’ed that the primary measure of the quality of the 
orbits derived from satellitc-to-satellitv Uueking would be obtained from orbit overlap test. 
Hence the orbit overlap test for the ATS-6/NIMBUS-6 experiment was performed in a way 
v/hich was more rigorous and less optimistic than tlic overlap test performed for tlie ATS-6/ 
GEOS-3 experiment. Notice that for the ATS-6/GEOS-3 experiment the overlap test was 
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ivrlormcil with data scit wliich intersected tlirouith the entire overlap interval, lienee both 
orbits used in the coinp.irison were constrained by data at each end ol the interval. With 
such a prircedure it is possible lor the el'lects of errors in the measuring system to cancel in 
the lest results. It will be* seen that lor the A I S-b/NIMHlIS-t) overlap test the two data sets 
•n question are abiding rather than overlaping and elfect ol measurement system errors are 
less likely to cancel in the lest results. 

Hie tracking schedules and the tracking systems used in l!ic evaluation are shown in 
t'igure ^). I he f irst two rows of this figure show the tracking schedules for the ranging data 
from Madrid, Spain to A l S-ti, and from Ahmed.ibad, India to ATS-ft. The third row 
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Figure 9. Meusurement Periods and Period for Nimbus Orbit Comparison 


indicates the tracking schedule for trilateration data (Madrid*ATS-6-Asccnsion Island). The 
fourth row shows the tracking schedule for the range sum and range sum rate data with 
Madrid as the ground station. The vertical bar located at 9HR/UT Feb. 8 indicates the 
e|'K)ch time for two estimated NIMBUS-6 epoch states. Kpoch state I was obtained by exe- 
cuting a Bayesian least scjuares estimator with all the ATS-6 ranging and trilateration data 
and all the satellite-to-satellite tracking data to the left of the epoch time. Epoch state 2 
was obtained by reflating the procedure with the satellite-to-satellite tracking data to the 
r«ht of the e| och time replacing the data to the left of the epoch time. The horizontal bar 
in row 5 of the figiu’e displays the common interval over which the two NIMBUS-6 epoch 
states were propagated. The complete GFM-7 gravity field model was used in all data reduc- 
tions and propagations. The relative weights for the data types weie obtained by first using 
nominal weights and processing all the data to estimate ATS-6 and NIMBUS-6 epoch states. 
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Tlic residuals of the cstimution were used to determine the "tandurd deviations of tlic noise 
oil the various data types. These standard deviations were used to obtain weights for the 
final data reductions. The computed standard deviations are shown on table 2. 


TABLE 2 

Standard Deviation Used for Wcigliting Measurements in Nimbus-6 
SatclIite-to-SatclIitc Tracking Orbit Determination 


RANGE (INDIA) TO ATS 

168 m 

RANGE (MADRID) TO ATS 

50 m 

TRILATERATION (MADRID, ATS, ASCENSION) 

15 m 

SATELLITE-TO-SATELLITE RANGE 

1 1 m 

SATELLITE-TO-SATELLITE RANGE RATE 

.3 cm/sec 


The data reductions were complicated by the fact that an experiment onboard the 
NlMBUS-6 was responsible for some outgasing. This effect was modeled as a constant along 
track thrust whose magnitude was estimated in the data reductions. The recovered magni- 
tude was approximately 1 0~7 m/scc2. 

Figures 10, 11, and 1 2 display the along track, cross track, and radial differences in 
the two epoch stale propagations during the overlaping period. The R.M.S. differences arc 
40 m along track, 30 m cross track, and 12 m radial. The secular growth of residuals in the 
along track direction is explainable in terms of gravity field error and an imperfect modeling 
of the outgasing effect whose direction was probably not exactly along track and whose 
magnitude was probably not constant. 

Finally it should be mentioned that a NIMBUS-6 orbit derived from satellitc-to-satel- 
lite tracking data was compared to a NIMBUS-6 orbit derived from minitrack data. The 
orbit differences were well within the stated accuracy for minitrack orbits of 500 m. 
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Figure 10. Along Track Differences for Nimbus-6 Satellite-to-Salellite 
Tracking Orbit Determination 

4.0 CONCLUSIONS 

The ATS-h/NIMBUST, and ATS-h/GfiOSO salclhUMo-aa.clli.e radio .racking sys.cm 

performs as specified wi.h a ,«olnlion of I meter in range and .03 em/sec in range-rale for 

a 1 0 second averaging. 

A Bayesian least sejuares estimation technujue utilizing a good a priori estimate of 
relay satellite state was used during these exf)criments to obtain user satellite orbits with 
accuracies comparable to what is obtainable from ground tracking systems. The limiting 
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Figure 1 1. Cross Track DilTerciices for Nimbus-6 Satellitc-U>-Sutcllite 
Tracking Orbit Determination 

factor in an orbit determination with satellite-to-satellite radio tracking appear- to be the 
accuracy of the force models rather than tracking system precision. 

The results of these experiments imply that with the proper data reduction proce- 
dures, the tracking data relay satellite system should provide orbit determination capability 
comparable to what is now obtainable from ground based systems. 
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